Extended Ginzburg-Landau formalism for two-band superconductors 
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Recent observation of unusual vortex patterns in MgB2 single crystals raised speculations about 
possible "type-1.5" superconductivity in two-band materials, mixing the properties of both type-I 
and type-II superconductors. However, the strict application of the standard two-band Ginzburg- 
Landau (GL) theory results in simply proportional order parameters of the two bands - and does not 
support the "type-1.5" behavior. Here we derive the extended GL formalism (accounting all terms of 
the next order over the small t — 1 — T/T^ parameter) for a two-band clean s-wave superconductor 
and show that the two condensates generally have different spatial scales, with difference disap- 
pearing only in the limit T ^ Tc- The extended version of the two-band GL formalism improves 
the validity of GL theory below Tc, and suggests revisiting of the earlier calculations based on the 
standard model. 

PACS numbers: 74.20.De, 74.20.Dw, 74.25. Ha 



The Ginzburg-Landau (GL) approach jjy, based on 
Landau's theory of second-order transitions, is one of the 
most powerful and most widely used theoretical tools of 
the present-day physics. It constitutes a solid base for 
theoretical studies in fields ranging from the condensed 
matter theory (e.g. superconductivity/superfluidity, 
phase transitions, and fluctuation phenomena) to par- 
ticle physics and cosmology (e.g. Higgs mechanism), 
and other topics reviewed in Ref. 0. It is generally be- 
lieved that the GL theory accurately describes essential 
physics in the vicinity of the critical temperature Tc (and, 
qualitatively, in a much wider temperature range). Sur- 
prisingly, this is not the case for two-band (and multi- 
band) superconductors, such as magnesium-diboride Q 
and several iron-pnictides Q etc., where the expected 
difference in spatial distributions of the two Cooper-pair 
condensates is not captured by the standard formulation 
of the two-band GL formalism. As explained further, the 
latter problem requires development of the extended GL 
theory, derived to a higher order in r = 1 — T/Tc- This 
is the core objective of this Letter. 

Recently, unconventional vortex patterns were ob- 
served in single-crystal MgB2 by Bitter decoration Q and 
by Scanning SQUID microscopy Although MgB2 is 
largely accepted as a type-II two-band superconductor, in 
Refs. ^ m no evidence of an Abrikosov lattice was found 
for low vortex densities. The interpretation was offered 
through the intervortex potentials derived from the stan- 
dard two-band GL theory of, e.g., Refs. Namely, 
for particularly chosen (different) coherence lengths S.i 
and penetration depths Xi of the two bands (i = 1,2), 
vortices were shown to conventionally repel each other 
only at short distances, while long-range attracting [lo| . 
This gives rise to stripe-like vortex patterns, unattain- 
able in either type-I or type-II superconductors, which 
led Moshchalkov et al. 'p\ to name this behavior "the 
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type-1.5 superconductivity". 

An avalanche of theoretical works followed 
based either on the standard two-band GL formalism 
itself, or the molecular dynamics simulations using the 
GL-calculated intervortex potentials, racing to describe 
the new type of superconductivity. Brandt was the first 
to point out that long-range vortex attraction is not nec- 
essarily a "type-1.5" property [l3|. The real criticism 
followed, in the analysis of Kogan and Schmalian fl^ . 
They showed that in the standard formulation of the 
two-band GL approach [i.e., two GL equations for two 
order parameters Ai(x) coupled through the Josephson 
intcrband coupling terms], there appear contributions 
to both order parameters of higher orders than r^/^, 
where t — 1 — T/Tc. However, the microscopic basis 
for the standard GL formalism assumes that only the 
terms oc r^/^ are accurate, which means that aforemen- 
tioned higher-order terms are incomplete and, thus, in- 
correct. After removing the higher-order contributions, 
Kogan and Schmalian found the order parameters of two 
bands to be proportional to each other and can be thus 
characterized by a single coherence length ^. As a conse- 
quence, type-1.5 superconductivity is not supported by 
this formalism. 

It is thus of abiding fundamental interest to clarify 
whether the relation Ai(x) oc A2(x) is generic to two- 
band superconductors or it holds only in the standard GL 
domain (to the order r^^^ in A^'s). To settle the above 
issues, we derive here the extended version of the GL 
formalism for a two-band clean s-wave superconductor, 
where the contributions to A^'s up to the order oc r^^^ 
are included in their full, correct form (while appearance 
of the higher orders is precluded). 

Our starting point is the BCS mean-field Hamiltonian 
of a two-band, s-wave, clean superconducting system. 
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i=i,2 



^],(x)t;-(x)^,,(x) 



+ ^]^(x)V;,Vx)A,(x)+h.c.J, (1) 

where j — 1,2 denotes each of the bands, H'^ is the c-term 
whose specific form (see, e.g., Ref. 0) is not of relevance 
for the present investigation, Tj(x) is the single-electron 
Hamiltonian, and the summation in the kinetic term is 
taken over the coinciding spin indices. The generalization 
of the mean-field self-consistency equation for two-band 
superconductors reads 



A,(x) 



E 

i=i,2 



(2) 



with gij being the relevant coupling constants [gij — gji). 
One of the most powerful formalisms to treat the su- 
perconducting properties in the presence of a nonuni- 
form spatial distribution of the pair condensate are the 
Gor'kov equations. For our study it is convenient to write 
these equations in the form of the Dyson equation for the 
2 X 2-matrix band propagator Qj^ (see, e.g., Ref. [iF 



where i?j(x) is a polynomial of Aj(x) and its spa- 
tial derivatives; = gijN{0) and nj = Nj{0)/N{0), 
where Nj{0) is the band-dependent density of states, 
and A''(0) = ^o construct the GL equations 

for a two-band superconductor, one should evaluate Ri 
with accuracy 0{t^^^). This results in two equations 
for Ai(x) and A2(x) coupled through the Josephson-like 
terms (for a clean two-band s-wave superconducting sys- 
tem, see, e.g., Refs.0,[3for a dirty two-band supercon- 
ductor, see, e.g., Refs.|8|, |9|). This is where the aforemen- 
tioned analysis of Kogan and Schmalian IJ] is important, 
as such a representation of the two-band GL equations 
must be corrected in order to avoid the appearance of 
terms of orders higher than r^/^ in A^ . Appropriate cor- 
recting procedure is given in detail in Ref. and results 
in two decoupled GL equations for Ai and A2 which ex- 
actly map on the one-band GL theory: Ai(x) cx A2(x) 
and both have the same coherence length unlike the ex- 
pectations based on the initial formulation of the two- 
band GL formalism. 

We now extend the GL formalism up to the order r^^^ 
in Aj's, by taking 
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where fiio = 7rT(2n -I- 1) is the fermionic Matsubara fre- 
quency (n is an integer and fcs is set to unity) and the 
2x2 matrix operator A^ in Eq. ([3]) is defined by 



Aj 
A* 



(x|Aj|x') =(5(x-x')Aj(x'). (5) 



Equations ^ and (0]) further give 

F- — C'-"'' A C 



!(0) 



q{0) a * ^^u; An, 



(6a) 
(6b) 



which makes it possible to expand J-j.uj in powers of Aj , 
when working near Tc- This is the well-known basis for 
Gor'kov's derivation of the GL theory 16 1. 

Using the definition of the anomalous (Gor'kov) 
Green's function 
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: — (X| J-,,, |X') = --(rV3,t(xt)V^,i(x't')), 



one can rewrite Eq. ^ in the form 

Ai(x) = Aiinii?i(x) + Ai2n2i?2(x), (7a) 
A2(x) = A2inii?i(x) + A22n2E2(x), (7b) 



with Af^ cx r^/^ and A^^^ oc r^/^. To begin with, we 
limit ourselves to a case of the zero-magnetic field (Aj's 
are real). Evaluating Rj with accuracy 0{t^^^), we ob- 
tain 



i?j = -aAj - bA^j -\- cA^ 

-f ICjV^Aj + QjV'^iV^Aj) - CjA.V ■ {Aj VA^), (9) 



with 



a = -M + r + -),yl = k(— --, 



b^Wi{l + 2T), w. 

93C(5) 



(1^3 - ^), 



(10) 



where hhjjj is the Debye energy, C,{. . .) is the Riemann 
zeta-function, F = 0.577 is the Euler constant, and the 
band-dependent Fermi velocity is denoted by Vj. Note 
that, as compared to the results of Refs. 7 and 1 14 there 
are three new terms in Eq. ([5]): oc A^, oc V^(V^Aj), and 

oc AjV • {Aj V Aj). In addition, the coefficients a, 6 and 
ICj contain extra contributions, i.e., a is now accurate 



up 



to the order whereas b and /Co include terms 



oc 



T. Note also that when evaluating Aj's with accuracy 
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©(r^/^) [see Eq. we have oc t in both Af^ and 

Going back to Eq. (O, one obtains Ri = Ai — 
Ai2?T.2-R2/(Aiini) from Eq. (I7al) . which can then be in- 
serted into Eq. (I7bp . Similarly, i?2 can be expressed as 
i?2 = A2 — A2inii?i/(A22'^2) from Eq. ((7b)) and substi- 
tuted in Eq. (j7ap . Such a manipulation, combined with 
Eq. results in the following equations 



aiAi 



a2A2 



biAl - ciAl - /CiV^A? - QiV2(V2Ai) 
+C1A1 V • (Ai VAi) - 7A2 = 0, 

b2Al - C2A\ ~ /C2V2A2 - Q2V2(V2A2) 

-h£2A2 V • (A2 VA2) - 7A1 = 0, 



where 



iV(0) 



-rynj(r + — ) 



r\ = A11A2 



^12' 



(11a) 
(lib) 

(12) 



with A\ — A22 — rjUiA and A2 = An — 77^2^1 {i] denotes 
the determinant of the A^ matrix, where A12 = A21). 
In addition, bj, cj, ICj, Qj, Cj in Eqs. (|llap and (jllbp 
are 6, c, /Cj, Qj, £j multiplied by njN{0), respectively. 
The last terms in the left-hand side of Eqs. (|llal) and 
()llb| are the Josephson interband coupling terms with 
7 = Ai2^(0)/r,. 

Proceeding in the manner similar to that of Ref. |l4 , 
we now group the terms of the same order in Eqs. (jllap 
and (jllbp . Keeping only terms of the order r^/^ in both 
equations we find 



aia2 
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= 0, 



(13) 



where (B)T-fe denotes the term in the expression B of the 
order r*^, with k an integer. Equation (|13p allows one to 
evaluate Tc in the two-band superconducting system and 
is reduced to A1A2 ~ ^12, which recovers Eq. (17) from 
Ref. H 

Further, when collecting the terms proportional to r"^/^ 
in Eqs. (fTTa]) and (jllbp we find 



aA 



(0) 



(0)i3 



KV^Af^ 



0, 



(14) 



where 





'<-{ 









/3i = 



6102 + a\b2/Y 



7 



/Cia2 -I- 1^20,1 
7 

/32 = /3iL^. 



, (15) 

where (52 is obtained from the expression for /3i by re- 
placing indices of a^'s and 6j's (1 o 2). Equation (jH)) 
is the correct formulation of the standard GL approach 
for the two-band s-wave clean superconducting system, 
as found in Ref. [l3- Using Eq. (fT3|). we indeed obtain 
from Eqs. dH]) and dH]) that 



[Ar(x)/Af(x)]^ 



A2/A1 



(16) 



which follows from the scaling /3i //32 = Ai/ A2- 

Now, taking the terms of order r^/^ in Eqs. (Illap and 
(jllbp . we arrive at 



Af\a + ?>P,[AfY)~KV^Af 



F{A!f)+F,{^"), (17) 



(0)n 



with 



F((^) =cr(^ + 5'VV + >'V^(VV), (18) 

and 

+ V^jVV^ + ^jfyJ^VV- (19) 

Equation (jl7p is the first main result of this paper. It 
includes all contributions to order t"^/^ to A^'s. Coeffi- 
cients (T, S and Y in Eq. (fT8|) are given by 



S 



Y 



aia2 
7 

Qia2 + Q2ai — /C1/C2 

7 



/Cia2 -I- /C2ai 

7 



(20) 



while the coefficients in Eq. (1191) read 

'6102 + af62/7^' 



Pi = 



Xi 



7 



C/i = - 



Cia2 - 30^6162/7^ -t- gf 02/7'^ 
7 

£102 + a? £2/7' 



61/C2 



7 

,^1-3 



a\lCib2 



(21) 



and /92, X2i C^2, and V2 are obtained from Eq. (I?T|) by 
replacing 1 •<->■ 2 in all relevant indices. 

Now, if the terms Fj{A^^^) were absent in Eq. (jT7)) . we 



.(1) 



would obtain that A^^'' (x) is proportional to A2^'' (x) and, 
furthermore, the ratio A^^''(x)/A2^''(x) would be identi- 



cal to Ai^^(x)/A^''^(x) given by Eq. (jT6| . However, in the 
presence of Fj{Af'^), this is no longer the case, as not all 



terms appearing in Fj ( A^"'' ) support the above scaling of 
the order parameters. In particular, let us consider the 
term /3j[A(*''']'^. This term could support the scaling only 
if the ratio pi/p2 is equal to Ai/A2- From Eq. (j2T|) we 
find 



pi _ Ai 2{niAl + n2Al) - rinin2{A2 + iAi] 



(22) 



P2 A2 2{niAl + n2Al) ~ i]nin2{Ai + 3A2) 
which means that pi/p2 7^ A1/A2 and, consequently, 

[AW(x)/A«(x)]V-42Mi. (23) 
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Moreover, as seen from the structure of Eq. p7|) . it is 



clear that A^^''(x) is not at all proportional to A2^''(x). 
We hereby arrive at our main conclusion, i.e., the band 
order parameters Ai(x) and A2(x) are not proportional 
to each other when extending the Ginzburg-Landau for- 
malism to terms in Aj's proportional to r'^/^ (beyond the 
standard terms oc r^/^). This means that the band coher- 
ence lengths are in general different, and this difference 
disappears only in the limit T ^ Tc- 

For completeness, we give here several remarks about a 
generalization of the extended two-band GL formalism to 
the case of a nonzero magnetic field (inclusion of a mag- 
netic field will not affect any of the above conclusions). 
Such a generalization is not straightforward because in 
the first step one needs to go beyond the eikonal approx- 
imation adopted by Gor'kov for the normal state Green's 
function (see, e.g., the textbook [17|). This task assumes 
extensive calculations with numerous details that are not 
suitable for a Letter. Therefore, we include here only the 
final result, while preserving the full derivation for a sep- 
arate publication: 



(x|dx') 



d 



a^ + i;-(---')i7^ ^(x|e-|x'), (24) 



' duj 



B^(x) 



where Gj^^ is the normal state Green's function in the 
presence of a magnetic field; the integration in the expo- 
nent is taken along a straight line connecting x and x'; 
(x — x')x is the component of the vector perpendicular 
to B(x) = rotA(x). As follows from Eq. (p4)) . the correc- 
tions to the Gor'kov approximation are gauge invariant 
and of order (A oc r^/^ and B oc r). In particular, 
Eq. ^ in the presence of a magnetic field reads 



3 h^n^^) 



-6A,|A,f + cAj|A,| 



^9 . ;x ^,^9^9 4m2r22(x) 



Aie 
3/ic 



rotB(x)D 



A - i 
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4|A,fD2A, 



3A*(DA,)2 + 2A,|DA,f + A^D^A, 



(25) 



where D = V -I- 27riA(x)/$o ($o is the superconducting 
flux quantum) and = |e|i3(x)/mc, with i?(x) = |B(x)|. 

As a final note, we state that our approach differs from 
the theory of a local superconductor in a slow varying 
magnetic field, used in Refs. 18, [l^ (the so-called gen- 
eralized Ginzburg-Landau-Gor'kov equations). The ap- 
proach developed in the latter papers assumes that the 
gradients of the order parameter are small but the order 
parameter itself can be close to its value at zero tem- 
perature. Instead, we extended the two-band Ginzburg- 
Landau formalism up to the order r'^/^ (in A^'s). This 



requires to accurately select the necessary terms on the 
basis of the proper scaling with r of A^ 's and their spatial 
derivatives. The same holds for the magnetic field and 
its spatial derivatives, which, contrary to Refs. 3 3 re- 
quires to go beyond the eikonal approximation of Gor'kov 
[see Eq. 

In summary, by developing the extended GL formal- 
ism for a two-band superconductor: (i) we improved the 
validity of the Ginzburg-Landau theory at temperatures 
away from Tc] (ii) we showed that the two position de- 
pendent order parameters in a two-band superconductor 
are generally not proportional to each other, thus their 
spatial scales are decoupled - contrary to conclusions of 
the standard GL formalism; (iii) we developed a useful 
tool for further theoretical studies of two-band supercon- 
ductivity, which also commands revisiting many earlier 
works based on the incomplete formulation of the two- 
band GL formalism. 
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